Dissipation, topology, and quantum phase transition in a one-dimensional Joesphson 
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We study the phase diagram and quantum critical properties of a resistively shunted Josephson 
junction array in one dimension from a strong couphng analysis. After mapping the dissipative 
quantum phase model to an effective sine-Gordon model we study the renormalization group flow 
and the phase diagram. We try to bridge the phase diagrams obtained from the weak and the 
strong coupling renormalization group calculations to extract a more comprehensive picture of the 
complete phase diagram. The relevance of our theory to experiments in nanowires is discussed. 
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I. INTRODUCTION 

Ohmic dissipation was shown to be of critical impor- 
tance in the macroscopic quantum coherence of a double- 
well system. ^'^ This zero-dimensional problem under- 
goes a spontaneous symmetry-breaking transition in the 
ground state if the resistor causing the dissipation is 
smaller than the quantum of resistance Rq = h/Ae^ ~ 
6.5kri; h is the Planck's constant and e the electronic 
charge. It is perhaps the simplest example of a quan- 
tum critical point (QCP).^ Soon afterwards, it was also 
shown* that a resistively shunted Josephson junction un- 
dergoes a similar phase transition from a metallic state 
(because there is always an Ohmic shunt present in the 
model) , in which the phase difference between the super- 
conductors is delocalized, to a superconducting state (in 
which a supercurrent can flow below the critical current), 
where the phase difference is localized in a minimum of 
the cosine potential of the effective classical Hamiltonian 
derived by Josephson. The root of this remarkable QCP 
is logarithmic interactions between the instantons^ rep- 
resenting the tunneling events, resulting in a dramatic 
failure of the dilute instanton gas approximation, so suc- 
cessful in the macroscopic quantum tunneling problemi^ 
The only difference between the two systems is the spe- 
cific ordering of the instantons;^ but in both cases the 
underlying cause is the orthogonality catastrophe® of an 
infrared divergent heat-bath that is necessary to model 
the Ohmic shunt resistor. Since then, numerous sugges- 
tions have been made that such an essentially dynamic 
zero temperature (T — 0) phase transition may be em- 
bedded in many condensed matter systems, and this may 
explain the ubiquity of dissipation in many ultra low tem- 
perature phenomena. Instead of enumerating here an in- 
complete list of references, we refer to Ref. ^ 

What about experiments regarding these predictions? 
We do not merely mean a reduction of the quantum 
tunneling rate due to dissipation, but a sharp phase 
transition or a QCP. Surprisingly, experimental evidence 
is sparse. The closest experiment that indicates sym- 
metry breaking in a double well is an indirect experi- 
ment involving a superconducting quantum interference 



devicei^ In contrast, in a more direct experiment on 
a single resistively shunted Josephson junction such a 
dynamic phase transition has been apparently experi- 
mentally observed. However, it seems that the the- 
oretical situation is more subtle than that assumed in 
the pastjM An attempt to observe such a transition in 
superconducting nanowires was also made*^ Unfortu- 
nately, the present experimental situation appears to be 
unclear iliSiiS-ii 

Sometimes the simplest theoretical concept is not the 
simplest from the experimental perspective. From the 
very beginning it was realized that such a dynamic tran- 
sition may have important consequences in many body 
problems^®' in particular in resistively shunted 
Josephson junction arrays (RSJJA). Beginning with the 
pioneering work of Orr et al?'^ experiments have been 
few and far in between. A brief recent review is pro- 
vided by Goldmani^ Of importance to us here are the 
experiments of Rimberg et al.r^ Takahide et al.r^ and 
Miyazaki et al?^ On the theoretical side, many impor- 
tant contributions have been made, but we list here only 
the pape rs that are germ ane to our present work; these 
are Refs. I27l28l29l30ll3lll3l33 ' 

RSJJA is a simple but nontrivial model, almost a 
paradigm to use a well-worn word, where both dynam- 
ics and statics can simultaneously play an important role 
in a quantum phase transition, as opposed to a classi- 
cal phase transition. Compared to a dissipative single 
junction problem, RSJJA is theoretically more challeng- 
ing because of the interplay of both spatial and temporal 
fluctuations. Indeed it is shocking that a recent theoret- 
ical work2M4 has found that the ground state of RSJJA 
is a state where the state of the (0 -I- l)-dimensional el- 
ements (single Josephson junctions) can slide past each 
other despite couplings between them. This is despite 
arbitrarily long-ranged spatial couplings. Such a phase, 
called the floating phase, was derived from a renormal- 
ization group analysis that is perturbative in the Joseph- 
son coupling. Given the striking nature of a lower- 
dimensional quantum criticality embedded in a higher 
dimensional manifold, it behooves us to examine the 
phenomena in the strong coupling limit and to see how 
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the two limits are reconciled. From the earlier hints, -'^^ 
the strong coupling analysis should show that when the 
Josephson coupling is gradually increased, the quantum 
phase transition changes its nature and cease to be en- 
tirely dynamic — topology, quantum mechanics, dissipa- 
tion, and the collective nature of the problem, all become 
equally important unlike the weak coupling limit. Mo- 
tivated by the theoretical challenge, recent experiments, 
and the paradigmatic nature of RSJJA, we analyze it in 
one dimension from the strong coupling limit. As the 
conceptually important issues are already present in the 
one-dimcnsional case, there is no need to examine more 
complicated higher dimensional situations. 

We use Villain mapping to investigate the strong cou- 
pling phase diagram and the critical properties. Villain 
mapping has been used previously to map the phase slip 
processes to a neutral gas of charges that have anisotropic 
logarithmic interactions in imaginary time and spatial 
directions because of dissipation i^^i^^i'^^'^'^W e improve 
upon the previous work and show that the resulting phase 
diagram is remarkably different. In a Josephson junc- 
tion array in the absence of dissipation, quantum phase 
slips are the topological excitations whose classical coun- 
terparts are vortices of the classical XY model. One of 
the basic concepts behind this mapping is a conservation 
law. Due to tunneling, the number of Cooper pairs will 
change on the superconducting grains. In the absence of 
any source or a sink this leads to a continuity equation 
on the lattice which serves as a constraint on the integer 
tunneling current. This constraint is resolved in terms of 
a single integer field defined on a dual space-time lattice, 
resulting in familiar vortices. Thus, the superconductor- 
insulator transition due to unbinding of vortices can be 
described by a one component sine-Gordon model. The 
presence of shunt resistances, or dissipation, alters this 
picture. Since we need to include the current through 
the shunts in addition to the tunneling current, we get a 
different conservation law with a source term (This point 
was missed in the past.) The presence of the source term 
does not allow us to resolve the current constraint in 
terms of a single integer field on the dual lattice, and we 
are left with a neutral gas of two-flavored charges and 
anisotropic long range interactions in space-time. Conse- 
quently, our phase diagram is different from those derived 
in the past. 

The rcnormalization group equations, derived here 
for the first time, are obtained by mapping the 
two-component neutral charge gas system to a two- 
component sine-Gordon model. A two component sine- 
Gordon model was used previously by Refael et al in the 
context of two resistively shunted Josephson junctions. 
There are therefore some similarities in the phase dia- 
grams. Apart from fully superconducting and metallic 
states some partially ordered phases are found. Surpris- 
ingly, the two-component sine-Gordon model also arises 
in the context of a classical two-dimensional XY model 
in the presence of an in-plane magnetic field.2i The mag- 
netic field spoils the zero divergence condition by intro- 



ducing a source term. Though it is clear that RSJJA and 
XY model in an in-plane magnetic field describe differ- 
ent physics, there are some similarities between the RG 
analyses of both problems. 

Our paper is organized as follows. In Sec. II we de- 
scribe the microscopic model and provide brieffy the weak 
coupling results. In Sec III we perform a detailed strong 
coupling analysis employing Villain mapping. To per- 
form RG calculations we map the two flavored neutral 
charge gas problem arising from Villain mapping to a 
two component sine-Gordon model. In Sec. IV we de- 
scribe the fixed point analysis of the RG recursion rela- 
tions. In Sec. V we construct the strong coupling phase 
diagram and contrast the results with those of the pre- 
vious authors. In Sec VI we briefly mention the relation 
between a superconducting nanowire and a RSJJA . Fi- 
nally we summarize our results in Sec VII. We present all 
the technical details of the RG calculations for the two 
component sine-Gordon model in the Appendix. Althogh 
this is relegated to the Appendix to preserve a smooth 
flow of the text, this Appendix is the heart of our theo- 
retical work. 



II. THE MODEL 

We consider the quantum action given by 

^ i ^ 'i,n 

- / drVy^cosid^e,) (1) 
Jo 

(3 is the inverse temparature, Eq = 2e^/C is the charg- 
ing energy, where C is the capaciatance of single island. 
ujn is the Matsubara frequency and a = h/Ae^R is the 
dissipation strength, where R is the shunt resistance. 
dxdi — Oi — ^i+i, and V is the Josephson coupling be- 
tween the grains. 

In the weak coupling regime <^ 1), which is dom- 
inated by strong quantum fluctuations, RG calculations 
lead to the following recursion relations upto order 0{V'^) 



dV_ 

Hi 
da 

ll 



{l--)V 
a 



^ = 



(2) 



So, a = 1 is a fixed line in the weak coupling limit. For 
a > 1, system is ordered (superconducting), and V is 
a relevant operator. For a < 1 the system is disordered 
(metallic), and V is irrelevant. Recently Tewari et al^^^ 
has extended this calculation to the order 0{V^) and 
have shown that for a one-dimensional array the above 
set of RG equations are correct upto third order. They 
also demonstrated that including longer ranged couplings 
that when V is irrelevant all longer range Josephson cou- 
plings are also irrelevant. They become simultaneously 
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relevant with V when a > 1. So, around the critical 
line a = 1, the chain of Josephson junctions decouples 
and behaves as independent single junctions. For these 
reasons the disordered phase (metallic) has been charac- 
terized as a floating phase. It is important to understand 
how far in the V ~ a plane is this weak coupling RG 
picture valid. 



III. A STRONG COUPLING ANALYSIS 

In the strong couphng hmit, that is V/Eq > 1, the 
phase difference between the neighboring grains will be 
localized around the minima of the periodic cosine po- 
tential and the tunneling events between the minima will 
govern the low energy physics. We analyze these tun- 
neling events using Villain mapping. To employ Villain 
mapping we first discretize the imaginary time into N 
slices so that f3 — NAt where At is the lattice spac- 
ing in imaginary time direction. After discretizing, the 
partition function becomes 

Z = J [DO] exp[J^ ^(1 - cosV ^9{l)) 



obtain 



-Y,^\u;\Arf{k)\0{q)\% (3) 



where the summation over repeated indices is implied, 
/i = (r, x) represents the components on the space-time 
lattice and I stands for the space-time coordinates of 
the lattice points. Here Jr — I/EqAt, = VAt and 
f{k) = 2[1 — cos(fca)]. For both Jr and Jx large we can 
replace the cosine by a periodic Gaussian, as follows: 

Z ^ J [DO] ^ exp[-:^X](^'^^«-2vrm^(r))' 



m^(i) I 
9 



(4) 



Here m^(i,r) are integers and represent tunneling be- 
tween minima of the cosine potential. 

Before we proceed further with the action described 
by Eq. ^ we explicitly demonstrate the difference of 
the conservation laws in non-dissipative and a dissipative 
Josephson junction array as have been mentioned in the 
Introduction. First we choose the non-dissipative case, 
and we set a = in Eq. 4. We also introduce two aux- 
iliary fields Bfj^ to decouple the periodic Gaussian terms 
as 



(V^0(/) - 27rm40))] (5) 



Summation over the integer field restricts the contin- 
uous fields Bf^ to integer values and integrating out 9 we 



Z cx ^(5v^B^,oexp[- 



vMl 



(6) 



So, we obtained the conservation law ^B^^ = in the 
form of a constraint on the configurations of B^. If we 
choose Bx = drh, Br — —dxh where ft is a single com- 
ponent integer field the constraint is resolved. For any 
nonzero a, we introduce an additional auxiliary field p to 
decouple the quadratic dissipative term: 



+ i 



(7) 



Now integrating out 9 we obtain the modified conserva- 
tion law 



V;.B^ + p = 



(8) 



which directly demonstrates that the zero divergence con- 
dition is destroyed and the auxiliary field p serves as a 
source term. 

We proceed with an analysis of Eq. 01 Integrating out 
9 from the partition function and doing the following 
transformations, 

Pi{i,T) = mr{i + l,T)-mr(i,T), (9) 
P2(«,t) = nixii.T + At) - mx{i,T), (10) 

we obtain 

Z = Zsw J2 exp[-2^' X]p^(g)G^,.(q>''(-<7)]. (11) 

In the above step we have rescaled time direction as 
At = AT^/VEoa. So, we have rewritten the parti- 
tion function in terms of two kinds of interacting charges 
Pi and p2 which by construction form a neutral gas of 
charges. Interactions of the charges are encoded in the 
Gfi^viq) which are given by 



Giiik,u; ) 

G22(fc,W ) 



J 



a (jj 



AT'a(fc2+w'2) 27rAr' (P -h w'2)' 

J a I 

AT'a(fc2 -t-w'2) 2^w'Ar' 

a 1 \ uj \ 



2tt At' (fc2 +u;'2)' 



J 



— G2I 



Ar'a(fc2 + w'2) 
In the long- wavelength and low- frequency limit, 
J „ „ a 1 



= (-1) 



i+j _ 



At' a{k^ + Lo'^) 



27r lu'At 



(13) 
(14) 

7, (15) 



4 



where J = J^J^ = jE^. From this interaction 
matrix it becomes clear that because of dissipation spa- 
tial kinks, i.e p2, have an additional anisotropic on-site 
logarithmic interaction in imaginary time in addition to 
isotropic logarithmic interaction in space-time. Charges 
of opposite sign for a given flavor attract and for differ- 
ent flavors charges of same sign attract. In the presence 
of dissipation we can not reduce the partition function 
in terms of a single set of charges or vortices given by 
Vi = V2 — V\ = £fiudfj,m^. So, we have to work with a 
two-component neutral charge gas, and the correspond- 
ing sine-Gordon theory will be a two component field 
theory. For convenience we shall choose Ar = 1/^/VEq, 
which leads to At = a. With this choice, the space- 
time lattice becomes isotropic and the anisotropy of the 
problem is captured only through the anisotropy of the 
interaction terms. 

We introduce two fugacities yi , y2 corresponding to 
Pi,P2 to control the charge fluctuations. As usual, for 
small 2/1, j/2 we can restrict pi,P2 to 0,±1. Thus, the 
sine-Gordon action is 

2 

k,uj i—1 

-^^y"dfcos(27r0i±27r02) (16) 

We are using the vector notation x = (x,r) and J dx = 
J dxdr. y+,y- are combinations of 2/1,2/2- At sec- 
ond order in fugacities, the RG procedures generate 
cos(27r(/)i ± 27r(/)2) terms, and that is the reason for ex- 
tending the coupling constant space, yi controls the for- 
mation dipoles between different signs of charges for i-th 
flavor. When the charges are bound as dipoles, the fugac- 
ities become irrelevant and flow to zero, implying order; 
when they become relevant free charges proliferate. 2/+ 
controls the density of dipoles between charges of same 
sign but of different flavors, and 2/- controls the density 
of dipoles between charges of opposite sign and different 
flavors; see. Fig. ^ is the matrix inverse of the in- 

teraction matrix for the charges. This 2x2 matrix is 
given by 



G7 



27r 



(17) 



Interestingly, the presence of divergent term in the 
propagator requires that we formulate the RG (thinning 
out the degrees of freedom) by integrating only in a given 
frequency shell. Then a dimensional analysis with the 
scaling prescription k — *■ ke~''^^,Lu — > cje"', where I is 
the length rescaling factor, shows that 



(18) 



as a is dimensionless. Since, |w| is nonanalytic and we 
have to maintain the periodicity of the cosine terms, the 




Quadrupoles 



> X 



FIG. 1: Configurations of space-time dipoles and qudrupoles 
of different flavors. Note that equivalent pictures can be 
drawn by charge conjugations. 



above RG equation for a is exact i^S*^ The following re- 
cursion relations up to second order are derived in Ap- 
pendix ^ 



dyi 
dl 

dy2 
dl 

dy+ 

dl 
dy_ 

dl 

dj 

m 



(1 + i - TT J)2/i + ay2y+ + (27r J + a)2/22/-(19) 
(1 + - - ^ J - a)2/2 + 27r J2/i2/_ , (20) 



-JT = (1 + a)y++'Kjyiy2, 



—— = (IH 27rJ- a)2/- - 7rJ2/i2/2, 



1, 



(21) 
(22) 



J[il--)-JiA,yt + A+yi + A^yi)]. (23) 



To keep the coefficient of fixed we need 



- = 1 - 

z 



(24) 



Here Aj^. A^ are regularization dependent constants. 
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IV. THE FIXED POINTS 



THE PHASE DIAGRAM 



Prom first order RG equations, tlie following picture 

emerges: j = l and the entire J — a plane is broken into 
six regions. In the region tt J > 2, a > 2 all the fugacities 
are irrelevant, and the system is fully superconducting 
and in the region (AttJ + a < 2) all the fugacities are 
relevant, and system is fully metallic. Other four regions 
are mixed phases where one or more fugacities become 
relevant thus implying special kinds of charge prolifera- 
tion processes. These four phases therefore have partial 
or mixed order. 

To find the fixed points of the second order equa- 
tions, we first note the structure of the equation for y+. 
Since, all the fugacities and J are always positive, when 
1 + - — a > 0, and Jj/ij/2 have to be zero. For Lorentz 
invariant fixed points corresponding to | = 1 and the 
coupling constant space specified byl + i — q;>0, which 
we can also write as A-^-y'^/2 + a < 2, we find the follow- 
ing fixed point solution: 

FPl: yi = y2 = 2/+ = y- = 0, J = J*, a = a* < 2. (25) 

When 1 -h i - a < i.e, ^+t/+/2 -h a > 2 , we get 

FP2: y^=y2=y+ = y_=0,J=r,a = a* > 2. (26) 

It is interesting to note that FP2 also corresponds to ^ = 
1. It is also worth emphasizing FPl and FP2 describe 
diff'erent parts of the J — a plane. They are not critical 
points; there are lines of critical points in the surfaces 
containing these fixed points, which are the same as those 
found from the first order RG equations. 

Consider now 1 + - — a = 0. For the results to be 

z 

sensible, z must be positive and therefore a > 1. The 
fixed points are now non-Lorentz invariant, and they are 
given by 

FP3: yi* = y*=y*_=0,r= 0, A+yf + 2a*=4 (27) 



FP4: yl = y^ =y*_=0, A+yf + 2a* = 4, J* = ^ 



(28) 



FP5: 2/2* = = 0, A+yf + 2a* = 4, J* = 



^2 (2-a*)(7r-2a*) ^ ^ tt 
'"^ M^* <2 ^^^^ 

These three sets of fixed points have continuously varying 
dynamic scaling exponents. However, despite the pres- 
ence of non-universal constants, z is universally deter- 
mined by a*. The only non- universality is in the loca- 
tion of the fixed points. Importantly FP4 and FP5 are 
the sought after intermediate coupling fixed points con- 
trolling the cross-over between the local criticality and 
the global criticality, which are not accessible in weak 
coupling calculations. 



In this section we construct the phase diagram, which 
is determined by FPl and FP2. This is illustrated in 
Fig. 2. In region, (tt J > 2, a > 2) all the fugacities are 
irrelevant and hence this phase is fully superconducting 
and is labeled as Superconducting. Here, each flavor of 
charges are bound in dipoles and the dipoles are formed 
between charges of opposite sign. Also, y+ is irrelevant, 
implying a quadrupolar order: dipoles formed between 
two different flavors of charges, but of the same sign, are 
controlled by y+ and hence in this region the inter-flavor 
dipoles are bound in quadrupoles. Following the simi- 
lar argument for we can conclude that inter-flavor 
dipoles formed between charges of different flavors and 
opposite signs are bound as quadrupoles. These different 
dipolar and quadrupolar orders are shown in Fig. 2. In the 
region, (47rJ + Q; < 2), all the fugacities, j/i , 1/2 , ?y+ , y- are 
relevant, and this is fully disordered or the metallic phase. 
Both flavors of charges proliferate and also the inter- 
flavor dipoles proliferate and we label this as Metallic. In 
the region, (ttJ > 2,a < 2), only yi,y2,y- are irrelevant 
and only y+ is relevant. This implies that each flavor of 
charges exhibit dipolar order but inter-flavor dipoles con- 
trolled by y+ proliferate. Hence, this is a mixed phase. 
This phase is labeled as Mixed 1. So, the transition be- 
tween the phases Mixed 1 and Superconducting is be- 
tween a dipole gas and quadrupolar order. In the region, 
(tt J < 2, a > 2) ,p2 charges are bound in dipoles and 
there is also quadrupolar order. But, pi charges pro- 
liferate. These facts are reflected by the irrelevance of 
y2,y+,y- and the relevance of t/i. For this reason this 
is a phase of mixed order and we label this as Mixed 2. 
In the region, (ttJ + a > 2,'kJ < 2,a < 2), y2,y- are 
irrelevant, and yi,y+ are relevant. p2 charges are bound 
in dipoles and inter-flavor dipoles controlled by y_ are 
bound as quadrupoles. But both the pi charges as well as 
the inter-flavor dipoles controlled by y+ proliferate, which 
is also an example of mixed order. This region is labeled 
as Mixed 3. In the region, AttJ + a > 2, ttJ + a < 2, 
only j/_ is irrelevant and all other fugacities are relevant. 
So, inter-fiavor dipoles controlled by y_ are bound as 
quadrupoles but pi,p2 charges and inter-fiavor dipoles 
controlled by y+ proliferate. This is also an example of 
mixed order and is labeled as Mixed 4. 



VI. CONCLUSIONS 

In this section we note our observations in regard to 
the recent weak coupling analysis of the fioating phase 
and briefiy describe the relevance of our theoretical work 
to various experiments that we have alluded to in the 
Introduction. In the near future, we hope to fully dis- 
cuss the experimental consequences of our theory. For 
nanowires, a number of important issues need to be ad- 
dressed, such as the effect of the external circuit for short 
wires, dissipative effects due to external circuit and the 
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FIG. 2: Phase diagram. Other than fully metallic and su- 
perconducting phases, there are four mixed phases character- 
ized by the fugacities that are irrelevant. The strong coupling 
phase diagram should not be valid for ttJ -C 1. The local crit- 
ical boundary at a = 1, at weak coupling, requires a separate 
analysis, as described in the text. 



resistance of the phase slip cores, boundary effects due 
to the leads^SiiMi and inhomogenitiesi^. For experi- 
ments involving arrays, the real time dynamics and the 
current-voltage characteristics would be a major theo- 
retical project. Presently, we can make only qualitative 
remarks that follow from our thermodynamic phase dia- 
gram. 



A. Floating phase and the local quantum criticality 

From the RG calculations at the second order, we 
found that 1/z depends on the fugacity y+. Therefore 
the locations of the nontrivial fixed points FP3, FP4 and 
FP5 found above depend on These fixed points can 
be accessed by tuning of in addition to tuning J 
and a. This is the reason why they do not appear in 
the phase diagram, which is a cut in the J — a plane. 
To recover the local quantum criticality observed in the 
weak coupling limit we note that as J is decreased be- 
low 2/7r, the fugacity yi becomes relevant, and we can 
no longer use the flows derived from small fugacity ex- 
pansions after a certain value of J. However, we observe 
that the growth of j/i implies that pi charges proliferate, 
and this fact allows us to restrict 0i « 0. Because of the 
anisotroipc imaginary time interaction, p2 charges may 



not proliferate at the same time. Using this approxima- 
tion, the two-component sine-Gordon problem reduces 
now to a one-component sine-Gordon problem described 
by the action 



k.uj 



— \lu\U2\^-2^ I dxcos{2Ti 
a J 



h) (30) 



Since, the propagator depends now only on Iwj, the 
problem is effectively zero dimensional. Note that this 
action by itself has an approximate duality with the 
weak coupling action in terms of the junction variables 
"01 = Qi ~ di+i under the identification y2 ^ V and 
a —> 1/a. The self duality at a = 1 is identical to the sin- 
gle junction problem. Thus, as the weak coupling limit is 
approached, at an intermediate coupling strength there 
is a change in the behavior from a global criticality to a 
local criticality. This can be easily recognized from the 
lack of self duality of the Josephson junction chain. By 
a simple RG calculation we will get 



dy2 
dl 
da 

m 



(1 - a)y2 




(31) 



If a > 1, 2/2 becomes irrelevant implying ordered state 
and as a < 1 the system becomes disordered. Here the 
ordering is a local ordering of the individual junctions. 



B. Josephson junction array 

In this subsection we describe the relevance of our the- 
oretical analysis regarding the experiments on resistively 
shunted Josephson junction arrays in one dimension. We 
have found that the phase boundary between supercon- 
ducting and fully metallic phase depends on the strength 
of the Josephson coupling of the superconducting grains. 
In our analysis we have predicted a complicated phase di- 
agram involving mixed phases in addition to fully super- 
conducting and fully metallic states, depending on the 
strength of J. Only two transitions (between Mixed 1 
and Superconducting and between Mixed 2 and Mixed 
3) take place which are independent of the strength of 
J. Experimentally determined phase diagram due to 
Miyazaki et al?^ confirms the strong coupling prediction 
that the phase boundary between ordered and disordered 
phase depends on J and it also shows that in the weak 
coupling limit there is a part of the phase boundary at 
a — \ which is independent of J . More detailed ex- 
periments are necessary to establish the complete phase 
diagram on a quantitative level for the entire J — a plane 
and would be possible when we establish the signature of 
many of the mixed phases, which are mainly dynamical. 

Apart from the region in which the system is fully 
metallic, all other regions correspond to either fully su- 
perconducting and mixed phases. We expect these super- 
conducting and mixed phases to demonstrate power law 
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behaviors for the temperature dependence of zero bias 
resistance: different mixed phases and the fuUy super- 
conducting phase can be distinguished by their different 
temperature exponents. Similar power law behaviors are 
also expected for the I-V characteristic. These power law 
exponents will depend on J and a which can be inferred 
from the linearized RG equations about the fixed points. 



C. Superconducting nanowires 

There is also a broad relevance of the dissipative 
Josephson junction array in the context of the experi- 
ments on superconducting nanowires pi2*i^4i^*iSiii AiT — 
quantum phase slips should play the key role in deter- 
mining if the system will be superconducting or resistive. 
An experiment by Bezryadin et al}'^ on a superconcon- 
ducting nanowire observed a dissipative phase transition 
similar to a single junction problem when the total nor- 
mal resistance of the wire exceeded Rq . Later worki^ 
implied that it is rather the resistance per unit length, 
hence the diameter of the wire, that is of significance. 

Quantum phase slips are topological excitations in 
which the phase of the superconducting order parame- 
ter slips by a quantized amount and can be viewed as 
a vortex in the space-time manifold at T = 0. If, for 
the moment, we can ignore dissipation and model the 
wire at T = as a (1 -I- l)-dimcnsional XY- model, it is 
well known that such a system cannot exhibit a sharp 
phase transition (Kosterlitz-Thouless) without consider- 
ation of topology or vortex unbinding; smooth Gaussian 
fluctuations, "spin waves" , cannot destroy superconduc- 
tivity. The finite temperature properties in the proximity 
of the QCP can be understood in terms of universal scal- 
ing functions, and for a finite wire, a finite size scaling 
analysis becomes necessary. 

For a quantum system, the above picture must be mod- 
ified because statics and dynamics are intricately inter- 
twined. So, one must specify an appropriate dynamical 
model. In one such model Cooper pairs are allowed to 
disappear in a pool of normal electrons.^" Thus, it is the 
actual phase of a superconducting grain that is coupled 
to an Ohmic heat-bath. While this may be sensible for 
a system of superconducting grains embedded in a nor- 
mal metal42i^, it does not seem to reflect the physical 
situation in a superconducting nanowire. 

The RSJJA is another model for a superconducting 
nanowire, at least as far as the global T = phase di- 
agram is concerned. The Josephson energy V provides 
an energy barrier for phase slips for currents up to the 
critical current, and the topological excitations that are 
so essential are automatically built in the model. At a 
coarse grained level, on the scale of a coherence length 
^, we can think of the wire to be partitioned into su- 
perconducting segments interrupted by the cores of the 
phase slips forming Josephson junctions. The character- 
istic frequency of the phase slips is given by the gap at 
which the conductivities of the normal metal and the su- 



perconductor are the same. So, the system will have three 
coupling constants a, Eq, and V corresponding to the de- 
gree of dissipation, the charging energy of the grains and 
the Josephson coupling strength. Here, the dissipation 
strength a — Rq/R^ and i?j = p^jixr^, where p is the 
resistivity and r is the wire radius. Importantly, we be- 
lieve that the correct model for the dissipative coupling 
is the one in which the phase difference across the phase 
slip is coupled to an Ohmic heat bath, as for RSJJA, 
unlike the model of Ref. fiol 

To determine the charging energy and the Joseph- 
son coupling strength, we consider the continuum action 
that describes the Mooij-Schon^ mode (gapless plasmon 
mode arising out of incomplete screening in one dimen- 
sion) is given by 

rL/2 . 
cfr / dx[cs{d^4>f + -{dr4>?] (32) 

2.TI Jq J-L/2 Cs 

where fi and Cs are the dimensionless kinetic admittance 
of the superconducting nanowire and the speed of propa- 
gation of the plasmon mode respectively. We have set 
h — 1. p and Cs are related to the superfluid den- 
sity and the capacitance of the wire by the relations 
^ E± ^ a.nd ^ ^ ^. A is the of the cross- 

27r 2 4m 27rCs 8e^ 

section, C is the capacitance per unit length, and the 
bulk superfluid density. If we discretize this action on the 
scale of a lattice spacing a, of the order of the supercon- 
ducting coherence length ^, we get the action of a Joseph- 
son junction chain with its parameters flxed by those 
of the Mooij-Schon mode. For the Josephson junction 
chain we get the following coupling constant relations: 
V = and ^ = For the strong coupling limit of 
the Josephson junction array the dimensionless param- 
eter J = fi/n. If we deflne kinetic inductance per unit 

length of the wire as i = 1/e'^ps we get p, — Rq/ \J L/C. 
Kinetic inductance per unit length can be expressed in 

terms of London penetration depth Al as L = Aig 
where Iq is the mean free path of the electrons. Con- 
verting this in terms of the normal resistance per unit 
length of the wire R we get L w i Tet^gT '^here we have 
restored h. 

We should note that, to leading order, we have J (x r, 
as the capacitance per unit length has only a weak loga- 
rithmic dependence on r. So, a phase transition depen- 
dent on J or p implies a dependence on the radius of the 
wire, that is, it tells us a critical radius beyond which su- 
perconductivity can be destroyed and these facts are in 
qualitative agreement with the results of Zaikin et alm^ 
and Biichler et al.'^^ . This fact is also in qualitative agree- 
ment with the experimental results when the resistance 
per unit length is plotted against the temperature>i^ In 
recent experiments on single-crystal Sn nanowires , Tian 
et al.^ has also emphasized the role of the diameter of 
the wire for very low temperature measurements, well 
below Tc- 

Upto now it might seem that there are no differences 
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between our results and those of Zaikin et al. and Biichler the quadratic order. If cjji has frequency and momen- 



et al. If we go through the data of Lau et al. care- turn components in the sheU a 



< uj < a 



and 



fully, we will find that for the MoGe wires used in the 
experiments a > 2. So, in these experiments only con- 
trol parameter is dimensionless admittance. This is the 
reason why our quantitative prediction fi — 2 is the 
same. But, significant differences will arise if the di- 
ameter and other system parameters can be adjusted 
such that a < 2; in this regime, we will observe tran- 
sitions between the mixed phases and also between the 
metallic and the mixed phases, depending on the cou- 
pling strengths. In the mixed phases, the resistance will 
follow various power laws as a function of temperature. 
For, a < 2, these exponents will be depend on both fi 
and a in contrast to the situation when a > 2, where the 
exponent depends only on fi. So, only if a as well as fi 
can be made sufficiently small, a local quantum critical- 
ity can be observed. Prediction of the new transitions is 
the outcome of treating the wire as an effective RSJJA. 
and we believe can be verified in more elaborate set of 
experiments on thinner wires and different materials. 

Before concluding this section we should mention that 
we did not take into account a few effects. One is the 
boundary effects due to the leads2S42i^ and the second is 
the possibility of inhomogeneties resulting in weak links. 
That inhomogeneties can play an important role is ev- 
ident from the recent experiments of Bollinger et alm^ 
Another experimentally relevant issue is the effect of dis- 
order on quantum phase slips. In a recent paper Khleb- 
nikov and Pryadkc""i considered the effect of disorder and 
demonstrated that disorder can bind the spatial coordi- 
nates of the phase slips and anti-phase slips and hence 
convert the problem to an effective (0-1-1) dimensional 
problem. They found that the phase transition takes 
place at fj, = 1 and belongs to the dissipative universal- 
ity class which we have defined to be the local quantum 
criticality. We hope to return to these interesting effects 
relevant to the experiments in a future work. 
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APPENDIX A: THE DERIVATION OF THE 
RENORMALIZATION GROUP EQUATIONS 



a/. 



< k < a' 



we define it to be the fast mode. 



otherwise, a mode is defined to be slow 



For a 



fixed frequency shell, we integrate out the the momenta, 
and then rescale both frequency and momentum in the 
resulting action for the slow part. The dynamic exponent 



to 



-jt/a ■ 



jt/a 



-(X)r 



FIG. 3; Renormalization group shell integration. The inte- 
gration proceeds by integrating out the shaded areas. 



z is necessary to capture the effects of the anisotropic in- 
teraction. The reason for this special order of integration 
is because of j^j in the propagator for ^2- If we reversed 
the order of integration, we would have encountered spu- 
rious infrared divergence from the frequency integration. 
While rescaling we must maintain the periodicity of the 
cosine terms, which implies that the field renormalization 
constant is unity. Cosine averages are being calculated 
with respect to the fast degrees of freedom being inte- 
grated out and it is denoted as ( )^ . For any operator 

O[(/ii,02], the average {O)^ is defined as 



{Oy= / D[$f]exp[~S[$f]]0[cf„,cl>2] (Al) 



In the small fugacity limit we can expand the par- 
tition function in powers of fugacities and truncate at 



We then arrive at the following action action for the slow 
degrees of freedom: 



9 



2 

/c,w i ■' r)=± ''^ j,j=l i=l,2,a=± a,/3=± 



where lij , lia and la/s represent the following second order contributions: 
kj = (^) / Jcos{2ncj>i{xi))cos{2wcj>j{x2W 



(A2) 



la,/} = ( 



(cos(27r0i(a;2) + 2'Ka(j)2{x2)) cos{2'!T(pi{xi)))-^ 



(cos(27r^i(fi) + 2'Ka(j)2{xi)) cos(27r^i(:r2) + 277/3^2(^2))) 



(A3) 



As previously used in the text, we are using the vec- 
tor notation x = {x, r) and = / dxdr. The calcu- 
lation involves the propagators G{j{x2 — x-i,T2 — ti) = 

G{j{x2 — xi) = G{j{r), which has an oscillatory behavior 
for long distance due to the sharp cutoff in frequency . 
With suitable regularization we can avoid the complica- 
tions that can arise duo to this oscillatory behavior. Dif- 
ferent regularizations lead to different nonuniversal con- 
stants in the recursion relations, but they do not affect 
the universal properties. We use the following regularized 
propagators 

G(,{r} = l-KoiX^/WT^) = -G{2(r-) = -G^^if) 

(A4) 



and 



GUr) = 5{x)£^Ko{X^f^^+^) + ^Ko{\yW+^) 



where A = o~^e~'. Kq is the modified Bessel function. 
For distance much greater than the cutoff these regular- 
ized propagators fall off exponentially and for distance 
much less than they have logarithmic dependence 
on distance. For example 



^ j dx{cos{2-iT(t>i{x)))f = ^ j dx{cos{2-K(t)l{x) + 2'k4{x))Y =^ J dxcos{2n^l{x))exp[-2TT'^{(f>'^i{x))f] 

= ^exp[-27r2Gf,(0)] J dxcos{27T^t{x)) (A5) 



After rescaling rescaling the coordinates we obtain 

yi{l) = exp[(l + -)l - 2n^GimyM (A6) 
z 

When we compute an average (cos(27r(^i(x)±27r(/)2(x)))-^ 
it involves computing the average {{tpi{x) ± (l)2{x))'^)^ . 
Due to the mutual interaction between and (^2 fields 
we obtain 

{mx)±4>2{x)fy 

= G{i(0) + G2^2(0)± 20(2(0) (A7) 

which leads to 

y±{l) = exp[(l-hi)Z-27r2(G{i(0)+G^2(0)±2G(2(0))]2/±(0) 

(A8) 



From equation (A5) and (A7) we obtain the first order 
RG equations for the fugacities. 



Now we illustrate the calculations of the second order 
terms. All the second order terms involve two space-time 
coordinates (xi,X2) but the propagators in real space- 
time are only functions of relative space-time coordinate 
{x2 — x\). For this reason we do a coordinate transforma- 
tion to the center of mass coordinate {R = 2i±i2.) and 
the relative coordinate {f = X2 — xi). The second order 
contribution In will involve only the isotropic propaga- 
tor G{j^(|r|). In order to calculate such a term we can go 
over to polar coordinates for the relative vector r which 



10 



has following three ranges 

df= d9[ rdr + / rdr + rdr] (A9) 

Jo Jo Ja Jae' 

In the first range of integration we can take the propaga- 
tor G(]^(|r|) = G(i{Q); in the second range either we can 
set G{]^(|r^) = G{i{a) or carry out the integral with the 
regularized logarithmic form of the propagator combined 
with a gradient expansion of the cosine terms. With suit- 
able regularization, G{j^(|?^) can be converted into a short 
range function of |?^ and third part of the integral can be 
made negligible The specific regularization controls 
how fast and how smooth the propagator falls off with 
the distance but does not affect the universal properties. 
This is what we achieve with the regularized propagator 
mentioned above. 

Due to the anisotropic part in the propagator G22 , the 
above decomposition into polar coordinates can be in- 

I 



First integral generates higher harmonic cos(47r0i(i?)) 
from the first integral range of the relative coordinate. 
The same integral in the second integral range generates 
a term proportional to the higher harmonic multiplied by 
square of the gradient of the field (j>i. So, obviously all 
these terms are irrelevant. Relevant contribution comes 
from the second integral when we concentrate on the sec- 
ond integral range of the relative coordinate. After a 

I 



As mentioned above here we have to deal with only an 
isotropic propagator 0(2 and hence polar coordinate will 
be useful. From the first range of the f integral we will 



correct for the second order terms involving the field (/>2 . 
So, we have to check all the second order terms men- 
tioned above independently. As we will demonstrate be- 
low the second order terms consisting of cross correlation 
of the fields , 02 lead to integrals over combination of 
G(i , G22 and 0(2 ■ If this combination turns out to be 
isotropic in space-time we can proceed with polar coor- 
dinate decomposition. Otherwise we need to follow a 
tricky decomposition which we will use below for /22. 
But, there is an interesting aspect of these second order 
contributions. If the second order term involves product 
of two different fugacities, it leads to relevant contribu- 
tions only if the combinations of the fields are contracted 
at the same space-time point which renormalizes other 
fugacities. This aspect simplifies the calculation of the 
second terms involving product of two different fugaci- 
ties. 

For illustrative purposes of the above procedures first 
we will pick In. 



(AlO) 



gradient expansion for cosine term we obtain, 

hi « -^yf J df(V0i)2 (All) 

Ai is a regularization dependent constant which in our 
case is Att^. This will contribute to the renormalization 
of J. Now we will pick a cross term e.g. 



(A12) 



get the contribution 

1,2 « (2^)e-2-'(G{,(0)+G4(0))(±4^2W (0)) X 

/ cos(27r(/)i(i?)T27r(/)i(i?)) (A13) 



2 /. 

= ^e-^"'°"("' / [(e~4"'^"(^"^ -l)cos(27r(0i(i? + r72) + <^i(i?-f/2))) 



R.f 



+ (e4-'G{,(r) _ ^) cos{2tt{MR + r/2) ~4>i{R- r/2)))] 

I 



1,2 = ^^^^ / /pi(27r0i(i?+r/2)+27r02(i?-r/2)) _^ ^ ^ j_ ^j(27r0i(i?+r/2)-27r02(i?-r/2)) 
20"* JR,? 



yi|2g-2.^(G{,(0)+G4(0)) f [(g-4.^G{,(f) _ 1) cos(2^0i r72) + 27r02 " r/2)) 



R,r 



+ (e4-'Gi2(P) _ ^) cos(2^,^i(i? + r/2) - 27r02(i? - r/2))] 

I 
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which renormahzcs y±. In the second range of the in- 
tegrals we will do a gradient expansion for the cosine 
terms which leads to the contributions proportional to 

cos(27r(/)i(^)±27r(/)2(^))x(VK(?!'iTVK(/'2)^- After Fourier 
transformation this leads to the terms which are combi- 



nations of frequency or momentum with cosine and hence 
irrelevant. 

Two interesting second order terms are /i+ and 7i_. 
For relevant part of these terms we get 



11+ 



h- 



1+1/^)1 f (e4^'(Gfi(rl+G{,(r)) _ cos(27r.^i (fl - f/2) + 2TrMR - ^/2) - 2nMR + 
Jr,? 



-2(1-1-1/0); 



(e4-'(G{i(fO-G(,(r)) _ cos(27r<^i(i? + f/2) + 27r02(i? - f/2) - 27r0i(i? - f/2)) 



(A14) 



Recalling that G(i + 0(2 = we can see that relevant /ii./i2,/i+ and exhaust all the second order terms 

part of I1+ vanishes identically. As both and that involve the fugacity yi. Now we will also exhaust all 

are isotropic in space time polar decomposition is valid the second order terms that involve the fugacity 2/2. For 

for Finally we get the relevant part of I2+ and h- we get 

/i_ w ^^47rJ« j dxcos{2'K(j)2) (A15) 



12+ 



y2{i)y+{i) 



-2{l+l/z)l 



J _ y2{l)V-[l) ^-2{\+\lz)l 
12— ~ A 



[ (e4^'(G4«+G{.(f)) _ 1^ cos(27r.^i(^ - f/2) + 27r</>2(^ - f/2) - 27r.^2(^ + f/2)) 

jR,f 

I (g4^^(G,^,(r)-G{,(r)) _ cOs(27r(/>2 (i? + f/2) + 27701 (i? " f/2) " 27702 (i? " f/2)) 
jR,r 

(A16) 

I 

Both of the combinations G22(^ ^ ^12 involve Finally we consider 722 which will exhaust all the terms 

f 

anisotropic part of the propagator G22- Since, the rel- 
evant contributions from the two terms above will come 
from the contraction of the fields at the same space-time 
point we can set f = without further difficulty and 
obtain 

I2+ « ^5^20/ j cos(27r0i) 
72- « ^^2(a + 27rJ)/ /dxcos(27r0i) (A17) 



that involve y2- 
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2 p 

^g-4^^G4(0) / [(e-4-'G4(r)_^)^Qg^27r(/.2(^ + f/2) +27702(7? -f/2)) 
0- jR,r 

+(e4'^'G4« _ ^-^ cos(27r02(^ + f/2) - 27r02(^ - f/2))] 

I 



(A18) 



Regardless of the anisotropic interaction the first integral Calculation of the second integral is tricky but leads to 
leads to higher harmonics of 02 and hence can be ignored. 
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contributions involving {dr4>2Y ^"^^ {dx4'2)'^ ■ But, in the 
original action such terms are not present and these are 
also of higher order than |w||02p and hence irrelevant. 
Though the second integral produces irrelevant terms we 



will describe the way to calculate it as such a calcula- 
tion will be required later on for / We will break 

the integral into two parts corresponding to only onsite 
contraction and ofFsite contraction of the field 4)2 ■ 



J 



122 



yl-A.-GUa) f ^^4.'GUr) _ 1) cos(2^</.2(i? + f/2) - 2^</.2 (i? - r/2)) 

a jR,f 

Ih^-in^aUo) / dxdudrie^'''^^^^'-^'^^ ~ 1) cos(27r</.2(x, u + t/2) - 2T:(f>2(x, " " r/2)) 
a-^ J 

^y|g-4.^G4(0) f ^^i^'Ggm _ ^) cos(2^02(i? + r/2) - 2^02(i? - r/2)) 

_y2g-4^^G4(0) / da;dudT(e'*'^'^2^(o,r) „ i)cos(27r02(2;,u + T/2) - 27r02(a;, u - t/2)) (A19) 

O'^ J 

I 



where 



u = '^^'Y^ and T = T2 — Ti. We have also bro- 
ken 6*22 into two parts 6*22 and corresponding to 
isotropic and anisotropic interactions respectively. Rear- 
rangement of /22 in this form has been used previously 
by Bobbert et aXwS For the first and third integrals which 
involve only onsite contraction we can break up the inte- 
grals into three parts 



dT 



dr 



dT 



dT 



(A20) 



Second range of the integrals contribute terms propor- 
tional to {du<p2)'^- It turns out that the contributions 
from the first and third integrals cancel each other. Sec- 
ond integral involves only isotropic function of |?^ and we 
can use polar coordinates as mentioned above which will 
contribute a term proportional to {V R(t>2)'^- 



Finally we have to consider only two other second order 
terms /++ and / Relevant part of these are given by 



^+(0 „-2(i+i/z); 



X / (e^"'(<=ii("^+^2'2«+2G(.«) _i)cog(2^((^i(i?-f f/2)-K02(i? + r/2)-</)i(i?-f/2) 



R,f 



X / (e4-'«=^i«+G4«-2G{.«)_i)cos(27r(0i(i? + r72) 

JR,r 



MR + r/2) -MR- r/2) + MR- r/2))) 

(A21) 



GUf)+Gi2{r) + 2GUr) 



if _ 5^x)^Ko[\^' 
depends only on time. This implies for /. 



^) which A+ is a regularization dependent constant and for the 



-i-i- 



the relevant 

contribution is obtained by contracting the combination 
of the fields 0i -I- (j)2 at same spatial point but at different 
times. So, we need to integrate over only t — T2 — ti. 
After breaking up the integral into three parts as we have 
done before the relevant contribution will come from the 
second range of the integral and we obtain 



-^Vl J dx{drM^ (A22) 



regularization used above, A+ = An . We have ignored 
the irrelevant terms like {9t02)^ and <pi(f>2 involving spa- 
tial or temporal derivatives. 

Since, G{^(r)-|-G22(r)-t-2G'{2(r) involves both isotropic 
and anisotropic parts, the evaluation of / is more com- 
plicated. But, we can carry out the calculation using the 
trick mentioned above for I22 ■ Ignoring irrelevant contri- 
butions we get 

/— « -^y- I dxi^M^ (A23) 
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where A- is another regularization dependent constant 
and for the regularization method chosen above A_ — 
47r3. 

Since, we are integrating out the momenta for a fixed 
frequency shell, the coefficient of / dx{dr4>i)'^ has to be 
kept fixed. Due to this reason after collecting the sec- 
ond order terms involving J dx{dr4>if' and rescaling the 
space-time coordinates we obtain, 

1 = + {A^yl + A+yl + A_yl)l] (A24) 

which leads to 

^ = - -) - J{Aiyl + A+yl + A^yl)] (A25) 



The anisotropic scaling prescription, k — > ke~^l^ can be 
used to keep the coefficient of fc^ fixed and 1/z is deter- 
mined by the difference of the coefficients of the terms 
involving J dx{dr<\>\f' and / dx[dx4>\Y generated at the 
second order RG transformation. This trick to extract 
the dynamic exponent in a perturbative analysis is well 
known. Following this trick we get 

2(1 - -) = A+yl (A26) 



This completes the derivation of the recursion relations 
of the coupling constants upto the second order. 
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